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ABSTRACT 


The  energy  spectrum  of  diffuse  neutrons  in  the  thermal  column  of  the  MIT 
reactor  has  been  determined  by  analyzing  transmission  data.  The  results  are 
compared  to  a  MaxMvellian  distribution  and  to  a  previous  experiment  done  on 
collimated  thermal  neutrons. 
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ABSORPTION  ANALYSIS  APPLIED  TO  NEUTRONS  IN  A  THERMAL  COLUMN 


INTRODUCTION 

It  has  been  shown  (see  reference  1  and  Appendix  1),  that  the  energy  spectrum 
of  a  collimated  neutron  beam  can  be  derived  from  transmission  data.  This  report 
will  describe  the  application  of  this  method  to  the  more  general  case  of  diffuse 
neutrons. 


THEORETICAL  PROCEDURE 

The  case  of  a  collimated  neutron  beam  will  first  be  discussed  and  then  the 
appropriate  changes  necessary  to  convert  this  to  the  diffuse  case  will  be  developed. 
In  reference  to  Figs.  1  and  4,  if  f  (  E )  represents  the  neutron  flux  per  energy  inter¬ 
val,  cr{E)  the  activation  cross  section  for  the  gold  detector, Sj  (E)  the  total  cross 
section  for  the  absorbing  material,  and  X  the  absorber  thickness,  then  the  activity, 
A(x)  ,  imparted  to  the  gold  detector  would  be 

CD 

A(x)=Kje"^T^^^V(E)f  {E)dE 

o 

where  K  is  a  constant. 

Figure  1  shows  the  geometry  used  for  the  diffuse  case.  Graphite  spheres 
borated  with  varying  amounts  of  boron  were  used.  Therefore,  instead  of  increasing 
the  thickness,  X  ,  of  the  absorber,  the  total  cross  section  2-p  increased  by 
adding  more  boron. 


The  Experimental  Setup 


As  can  be  seen  from  Fig.  1,  a  neutron  may  reach  the  gold  detector  by  more 

-S  (E)x 

than  one  path.  Hence,  the  attenuation  factor  which  was  6  T  for  the  beam 

case  is  now  more  complicated,  namely, 

r-5:T(E)x 

e  '  P(x)dx  (2) 

I 


where  P  (  X  )  is  the  fraction  of  neutrons  that  travel  the  distance  X  per  X  interval. 
Fortunately,  it  was  found  (see  Appendix  2)  that  this  attenuation  factor  could  be  repre¬ 
sented  to  a  good  approximation  by  a  simple  exponential, 
oo 


-2t(E)x  -2t(E)>< 

e  >  P(x)dx  =  Ke  ^ 


(3) 


I 

where  "x  can  be  considered  as  an  equivalent  distance  through  which  the  average 
neutron  passes,  and  K  is  a  constant. 

Now  Eq.  (1)  can  be  rewritten  for  the  diffuse  case  as 


where 


oo 


A(B)  =  K 


e‘^®cr(E)f(E)dE 


2  = 


0-0  (E)N 


(4) 


(5) 


B  =  ^  X  b  (6) 

where  the  quantities  O'0(E),N  ,  A0  ,  and  b  are  defined  as  the  microscopic 
total  cross  section  for  boron,  Avogadro's  number,  the  atomic  weight  of  boron,  the 
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sphere's  density,  and  the  weight  fraction  of  boron,  respectively.  It  should  be  noted 
that  the  quantity  B  has  the  dimensions  of  grams  cni~  ,  which  makes  it  suitable  for 
a  transmission  curve,  and  also  tliat  X  does  not  h:ive  its  usual  meaning  and  is  inde¬ 
pendent  of  the  boron  concentration. 

Now,  since  X,  ®’(E),  and  f  (E)  are  all  functions  of  E  ,  we  can  define  a  new 
quantity,  •pa)  ,  such  that 


^  (2)d2  =  <^(E)t(E)dE  (!) 

If  we  now  substitute  Eq.  (7)  into  Eq.  (4),  we  obtain 

ao 

A(0)  =  -K|e'^®<#>(S)di:  I®) 

o 

which  involves  the  assumption  that  X  -  0  for  E  ,  and  X  =  oo  for  E  *  0. 

Equation  (8)  now  has  the  form  of  a  Laplace  transform  and  hence,  if  A(B)  is 
determined  experimentally,  can  be  obtained  from  the  inverse  transform  of 

A(B). 

The  desired  quantity,  f  (  E  )  ,  is  then  obtained  by  use  of  Eq.  (7); 


f{E)  = 


<^(E) 


(») 


EXPERIMENTAL  PROCEDURE 

Figure  1  shows  the  experimental  setup.  A  gold  detector  was  placed  in  the  center 
of  the  two  borated  hemispheres.  This  was  held  in  a  graphite  cylinder  which  was  inserted 
into  a  graphite  stringer  for  the  irradiation  with  thermal  neutrons  at  the  MIT  reactor. 

Six  inches  above  the  sphere,  another  graphite  cylinder  containing  a  gold  detector  was 
inserted.  This  acted  as  a  flux  monitor. 


-  4  - 


The  activity,  A(B)  j  of  the  gold  ball  within  the  i*’^  borated  sphere 
would  be  given  by 


00 

A(B)i  =  N(B)jSj  §{B). 

o 

where  N  (B)j  is  the  number  of  gold  atoms  in  the  i^^  detector,  Sj  is  the  activation 
saturation  term,  and  $(B)j  is  the  flux  at  the  detector. 

The  activity  A(m).  of  the  i^^  monitor  would  be 


(10) 


A(m)j=  N(m)j  Sj  $(m)j  o- 


(11) 


where  N  (m)  |  is  the  number  of  gold  monitor  atoms,  $  (m)j  is  the  flux  at  the 

monitor,  and  a  is  the  average  activation  cross  section  for  gold. 

What  we  actually  measured  was  the  ratio  of  Eq.  (10)  to  £q.  (11)  which  would  be 

CO 


A(B)i 

A(m)j 


=  R{B)  = 


N(B)jSi$(B)j  Je 

N(m)  I  Si  $  (m)j 


(12) 


Equation  (12)  simplifies  to 


R(B)  = 


W(B)i 

W(m)| 


4>mdl, 


(13) 


W(B)|  ® 

where  - -  is  the  weight  ratio  which  would  be  equal  to  the  atom  ratio,  and 

W  (m)j 

F  is  the  ratio  of  the  fluxes  with  no  boron  present 


$(B)o  A(B)o  W(m)o 
f(m)o'w(B)o  A(m)o 


Equation  (1^  can  be  rewritten 


R(B)W(m)i 
F  W(B)i 


oo 


=  P(B)  = 


(14) 


(15) 
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With  no  boron  present,  P(B)  ^vould  be 


P(0)  = 


R  (O)W(m  )o  a  _ 
F  W(B)o 


Therefore,  in  order  to  make  P(0)=  I  ,  we  define  Q  (  B)  as 


where 


Q(B)  = 


P(B) 


oo 

r 

-S  B 

e  0(5:)d2 

o 


S(Z)  = 


Finally,  the  desired  quantity, 


<^(Z) 

a- 

f  (E)  ,  would  be  given  by 


<16) 


(17) 


(18) 


f(E) 


9{^)cr  d  Z 

^^(E)  Te 


(19) 


where  eil)  is  the  inverse  Laplace  transform  of  the  experimentally  determined 

quantity  Q(B),  and  the  other  factors  are  obtained  from  known  cross-section  data  for 
2 

boron  and  gold. 

As  can  be  seen  from  Eq.  (17)  and  (15),  one  needs  to  know  the  ratio  of  the  activi¬ 
ties,  R{B)  ,  in  order  to  calculate  Q(B).  This  was  determined  by  placing  the 
irradiated  detectors  and  monitors  into  a  scintillation  well  counter  and  measuring  the 
ratio  of  their  respective  count  rates,  which  should  equal  the  ratio  of  their  activities. 
This  ratio  should  also  be  independent  of  time. 


RESULTS 


Table  1  shows  Q(B)  along  with  the  values  of  B  for  the  various  spheres  (see 
Appendix  3).  The  transform  T(B)  yielded  the  best  fit  to  Q  ( B )  •  T( B  )  is  given  by 


T(B) 


-  Jj) 


(20) 
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where  0  and  d  were  found  to  be  16.  76  and  0.  0256,  respectively.  This  transform 
is  also  found  in  Table  1. 


The  inverse  transform  of  Eq.  (20)  is 

Q\f6 


ae 


2i/Tf 


-dS  -3/2 

e  e  .  T 


(21) 


Now  f  (  E)  can  be  obtained  from  Eq.  (21)  with  the  use  of  Eq.  (19).  This  is  shown 
compared  to  a  Maxwellian  flux  distribution  at  293°K  in  Fig.  2. 


CONCLUSION 

Previous  experience^  has  shown  that  this  method  should  be  valid  within  10 
percent.  However,  as  can  be  seen  from  Fig.  2,  the  experimentally  determined 
spectrum  differs  from  the  Maxwellian  distribution  by  greater  than  10  percent.  There 
are  at  least  four  possible  explanations  for  this  fact; 

1.  The  neutrons  in  the  MITR  thermal  column  have  the  energy  distribution  as 
found,  within  10  percent. 

2.  The  mathematical  analysis  of  the  problem,  especially  the  scattering  part 
thereof,  may  be  in  error  to  a  degree. 

3.  The  positioning  of  the  spheres  relative  to  the  monitors  was  not  done  with 
sufficient  accuracy. 

4.  The  flux  gradient  between  the  sphere  and  the  monitor  did  not  remain 
constant  throughout  the  experiment. 

Of  the  possibilities,  3  and  4  seem  the  most  probable,  since  this  experiment 
was  carried  out  over  a  three-week  period  which  involved  ten  separate  irradiations. 

However,  it  appears  that  the  method  of  absorption  analysis  can  be  applied  to 
the  diffuse  case,  and  it  certainly  could  be  improved  by  changing  the  experimental 
setup  so  as  to  ensure  constant  geometry  and  perhaps  eliminate  the  scattering 
problem. 


-  7  - 


VELOCITY  IN  METERS  PER  SECOND 
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APPENDIX  1 


ABSORPTION  ANALYSIS  ON  A  COLUMATED  NEUTRON  BEAM 


For  the  collimated  case,  silver  absorbers  of  thickness  X  ^vere  used,  with  a 
boron  trifluoride  tube  as  the  neutron  detector.  The  absorption  curve  obtained  is 
shown  in  Table  2  along  with  the  transform  fit.  The  transform  that  yielded  this  fit 
was 


1.1 


Equation  1. 1  yielded  the  inverse  transform, 

0  when  (2-a)<0 

2 

b  ^  [byc-c(S-o)-^^|:^]  1.2 

2Tr(5:-a)3/2  »heD  (S-a)>0 

where  the  values  of  Q  ,  b  and  c  were  found  to  be  1. 30,  3. 31  and  0.  31,  respectively. 

In  order  to  obtain  the  neutron  flux  distribution  from  4>{1  ),  it  must  be  multi- 
ri  5* 

plied  by  ■  H  and  divided  by  the  energy-dependent  detector  efficiency.  The  result 

d  V 

thus  obtained  is  shown  compared  to  a  Maxwellian  at  293*^K  in  Fig.  3. 


0(2)= 
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VELOCITY  IN  METERS  PER  SECOND 


APPENDIX  2 


DETERMINATION  OF  THE  SCATTERING  CORRECTION  FACTORS 


Calculation  ^viU  be  made  of  the  effective  length  that  the  neutrons  must  pass 
through  to  reach  the  gold  detector  at  the  center  of  the  sphere.  The  results  obtained 
here  will  show  that  the  neutron  attenuation  is  exponential  with  respect  to  the  total 
cross  section.  Thus,  the  neutrons  act  as  if  they  were  travelling  through  1.062  cm 
oi  absorber  in  a  1. 000-cm  radius  sphere.  This  simple  exponential  property  allows 
us  to  use  the  Laplace  transform  technique  to  solve  the  integral  equation. 

Figure  4  shows  the  geometry  for  this  problem.  Let  us  consider  an  elemental 
volume  in  this  figure.  The  number  of  neutrons  that  reach  it  from  a  unit  area  on  the 
surface  would  be 


#(x,e) 


2.1 


where  ^  is  the  diffuse  flux  at  the  surface,  and  S  y  is  total  macroscopic 
cross  section  of  the  sphere.  The  quantity  ^  (  X  ,  ^ )  would  be  the  neutron  current  at 
the  point  X  for  a  given  0,  Now  if  Jg  is  the  macroscopic  scattering  cross  sec¬ 
tion,  then  <f>  {  X  , ^  )  would  be  the  number  of  scattered  neutrons  per  second  per 
unit  volume.’  If  we  now  multiply  this  quantity  by  the  number  of  unit  volumes  dv, 
we  have 


df  =  <#)(x,0)2s«iv 

where  d  f  is  the  differential  number  of  neutrons  being  scattered  per  second  per 
differential  volume.  Combining  £q.  2. 1  and  2. 2  we  get 

-2tX 

/  6-  \ 

df 


2.2 


-  ^ 
\  ^  /  2ttx 


211x2 


2.3 
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Figure  4.  The  Scattering  Geometry 


From  Fig.  4,  dv  would  be 


dv  =  x^sin^dx  d0d»// 

Combining  Eqs.  2.3  and  2.4,  we  get 


2.4 


“StX 


df  =  e  ^  2c  sin  0dx  d^  di// 
Qyr  3 


2.5 


The  neutrons  at  df  must  travel  the  distance  y  in  order  to  reach  the  center. 
Therefore,  the  differential  flux,  d  n  >  ihe  center  would  be 


n  rtf 

dn  = - ^ —  df 


4Try2 


From  Fig.  4 


X  ^+-  1-2  X  COS  0 


2.6 


2.7 


Now  by  combining  Eqs.  2, 5,  2. 6  and  2. 7  we  get 


dn  = 


-2-(  x+y  X ^  +1-2  X  COS0  ) 

T  c5 


sin^dxd^d^/  2.8 


32tt2(x^+1-2  XCOS0) 


In  order  to  obtain  the  total  flux,  ^  ,  at  the  center,  due  to  the  entire  surface, 
we  must  multiply  by  417  r  ^  =  4 17  and  integrate,  thereby  obtaining 


2.9 


2cos0  WZ  2TT 


N  = 


.  <^0^3 


817 


-Z.-(x+\/  x^-f- 1 -2  XCOS0)  . 

»  T  SI 


sin^  dx  d^  dv^ 


x=0  0=0  -0 


(  X  ^  +  i  -  2  X  cos  9  ) 
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Since  the  integration  with  respect  to  '/z  amounts  to  multiplication  by  2  7T  and 
Sg  =0.48  cm  ^  ,  we  now  have 


2cos8 


N=O.I2<#)o 


/  X  ^4- 1-2  X COS0  ) 


sin 0  dx 


2.10 


{  x^4- 1  -2  X  cos  9  ) 

x=0  0=0 

Since  this  N  represents  the  neutron  current  at  the  center,  for  activation 
purposes,  2  N  niust  be  used. 

Unfortunately,  Eq.  2. 10  cannot  be  evaluated  analytically.  Graphical  methods 
must  be  used.  The  results  of  this  double  graphical  integration  are  shown  in  Table  3. 
The  flux  at  the  center  due  to  primary  scattering  is  represented  by  repre¬ 

sents  the  straight  in  component,  and  the  flux  at  the  center  due  to  <f)^  and  <f>  ^  . 
The  value  of  at  (0.50)  corresponds  to  no  boron.  This  shows  that  the 
dlrect-in  plus  the  primary  scattered  neutrons  account  for  92.  3  percent  of  the  neutrons 
at  the  center.  The  remaining  7, 7  percent  would  be  due  to  multiply -scattered  neutrons. 
This  7. 7  percent  is  too  large  to  ignore,  and  since  it  cannot  be  calculated  conveniently 
at  other  values  of  Z  y  ,  it  is  estimated  to  drop  off  at  the  same  rate  as  the  primary. 
This  is  equivalent  to  multiplying  the  primary  contribution  by  1.  247  to  give  now  the 

total  contribution  at  the  center  due  to  scattering,  designated  as  tabulated 

bl 

along  with  j  j  which  is  now  +  <f>^ 

It  was  found  that  ^yy  could  be  represented  by  a  simple  exponential 


^^.^=*-'■062(2^-0.5) 


2.11 
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APPENDIX  3 


DETERMINATION  OF  BORON 

The  berated  spheres  were  made  by  mixing  known  amounts  of  reactor  grade 
graphite  powder  with  purified  boron  powder.  After  shaking  on  a  paint  mixer  for  15 
minutes,  the  samples  were  subjected  to  a  pressure  of  12  tons  per  square  inch  to  form 
a  solid  berated  graphite  rod.  From  this  rod  two  1-cm  hemispheres  were  machined 
and  used  to  hold  the  gold  ball  as  shown  in  Fig.  1. 

The  percentage  of  boron  was  determined  by  three  different  methods;  1)  the 
weights  used;  2)  chemical  analysis  of  the  mixed  powder;  and  ^  neutron  transmission 
analysis  of  the  hemispheres. 

The  chemical  analysis  was  made  by  converting  the  boron  to  boric  acid,  com- 
plexing  the  boric  acid  so  as  to  make  it  a  strong  acid,  and  then  titrating  with  standard 
hydroxide. 

The  neutron  transmission  analysis  was  carried  out  at  the  MIT  reactor  using  a 
collimated  beam  of  monoenergetic  neutrons  and  determining  the  percent  transmission 
for  the  various  hemispheres.  Since  these  neutrons  had  an  energy  of  0.0806  electron 
volts,  for  which  boron  has  a  total  cross  of  422  barns,  the  quantity  B  is  given  by 

B=  In —  =  0.04255  X  In  ~ 

No-  T  T 

where  A  is  the  atomic  weight  of  boron,  N  is  Avogadro's  number,  cr  is  the 
microscopic  absorption  cross  section  of  boron,  and  T  is  the  transmission. 

The  results  of  the  three  methods  are  shown  in  Table  4.  As  can  be  seen,  the 
claim  of  99.2  percent  for  the  purity  of  the  boron  seems  to  be  a  bit  high.  The  trans¬ 
mission  analysis  data,  the  only  method  which  actually  would  take  into  account  the 
cross  sections  of  any  possible  impurities  in  the  hemispheres,  was  used. 


TABLE  1 

VALUES  OF  Q(B),  T(B),  and  B 


Sphere 

B  (grams/cm^) 

Q(B) 

T(B) 

1 

0 

1.000 

1.000 

2 

0.00417 

0.  804 

0.814 

3 

0.0118 

0.579 

0.  574 

4 

0.0174 

0.448 

0.455 

5 

0.03099 

0.272 

0.271 

6 

0.0471 

0.  172 

0. 161 

7 

0.0655 

0.0776 

0.093 

8 

0.08232 

0.0593 

0.0600 

9 

0.1049 

0.0346 

0.0346 

10 

0.1296 

0.0205 

0.0199 
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TABLE  2 

COMPARISON  OF  T(x)  WITH  THE  EXPERIMENTAL  ABSORPTION  DATA 

FOR  THE  COLUMATED  CASE 


x(cm) 


0.0000 
0.0230 
0.0591 
0. 1322 
0. 1811 
.238 
.324 
.430 
.622 
.  795 
.952 
.  190 
.365 


Absorption  Curve 

T(x) 

1,000 

1.000 

0.  908 

0.  907 

0.783 

0.  783 

0.586 

0.586 

0.491 

0.491 

0.399 

0.  399 

0.297 

0.  297 

0.209 

0,209 

0.1152 

0. 1150 

0.0691 

0.0693 

0.0446 

0.0445 

0.0230 

0.0233 

0.0148 

0. 0148 

TABLE  3 


VALUES  OF  ^TT’^ST’ 


2  j(cm'‘  ) 

(/>S 

4>J 

‘^ST 

'^TT 

0.50 

0.316 

0.606 

0,923 

0.  394 

1.000 

1.20 

0,1344 

0.301 

0.435 

0.  168 

0.419 

2.50 

0.0275 

0.0820 

0. 1095 

0.0343 

0. 1163 

4,00 

0.00502 

0.0183 

0.0233 

0.00625 

0.02455 

6.00 

0,000556 

0. 00249 

0.00305 

0,000695 

0.00318 

TABLE  4 

VALUES  OF  b  FROM  DIFFERENT  ANALYSES 


Sphere 

Weights 

(Assuming 

0.  992  Purity) 

Chemical 

Analysis 

Neutron 

Transmission 

Analysis 

2 

0.00212 

0.0019 

0.00185 

3 

0.00556 

0.0049 

0.0052 

4 

0.00820 

0.0076 

0.0075 

5 

0.0151 

0.0144 

0.01374 

6 

0.0227 

0.0216 

0.02084 

7 

0.0322 

0.0306 

0.02945 

8 

0.0402 

0.0378 

0.03735 

9 

0.05145 

0.0505 

0.0483 

10 

0.0628 

0.0592 

0.0599 

-  19  - 
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